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We discuss experimental and observational constraints on certain models of higher dimensional 
light Majorana neutrinos. Models with flavor blind brane-bulk couplings plus three or four fla- 
vor diagonal light Majorana neutrinos on the brane, with subsequent mixing induced solely by the 
Kaluza-Klein tower of states, are found to be excluded by data on the oscillations of solar, atmo- 
spheric and reactor neutrinos, taken together with the WMAP upper bound on the sum of neutrino 
masses. Extra dimensions, if relevant to neutrino mixing, need to discriminate between neutrino 
flavors. 
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Do neutrinos have anything to do with large extra di- 
mensions Do these extra dimensions discriminate 
between neutrino flavors? These questions have engaged 
the attention of many authors [2| during the last few 
years. Let us consider the simple picture of a three-branc, 
containing the Standard Model (SM) fields, embedded in 
a higher dimensional space, known as the bulk. The ba- 
sic idea then is to introduce one (or more) right chiral 
neutrino(s), which are singlet(s) with respect to the SM 
gauge group, in the bulk [3j. Tiny physical masses for 
the active SU(2)l doublet neutrinos on the brane, that 
mix among themselves as well as with sterile components, 
then get induced. These are controlled by the size(s) of 
the compactified extra dimension(s) and the strength(s) 
of the couplings(s) between the brane and the bulk states. 
The bulk states expand into a Kaluza-Klein tower of 
states with a zero mode, as well as accompanying KK ex- 
citations of heavier sterile neutrinos with arithmetically 
progressing masses = ±k/R. k=l,2,... with R being 
the radius of the largest extra dimension. The leakage of 
a propagating active neutrino into these sterile states has 
to be kept controllably small so as not to significantly al- 
ter the conventional oscillation probability formula [5j in 
the standard picture of three directly mixed active neutri- 
nos. An immediate concern, nonetheless, is whether this 
picture is compatible with the pattern of neutrino masses 
and mixing angles implied by the latest oscillation and 
cosmological data. This note will address precisely this 
issue, focussing on the question whether the brane-bulk 
couplings are flavor blind or whether they are forced by 
all the data to have significant flavor dependence. 

A relevant basic issue is whether the physical massive 
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neutrinos are Dirac or Majorana particles. There is a rel- 
atively straightforward picture in the former case. One 
combination of the bulk neutrinos now becomes the right 
chiral partner of an active neutrino on the brane. To- 
gether they acquire a Dirac mass mp, ~ hvMp/Mpi, 
where v is the electroweak vacuum expectation value 
(VEV), h is a Yukawa coupling strength (controlling the 
interaction between the brane and bulk states) which 
is of order unity or less, Mpi the reduced Planck mass 
in (1+3) dimensions whereas Mp is the higher dimen- 
sional fundamental scale, proposedly in the multi-TeV 
range. This scenario was given a thorough phenomeno- 
logical examination in Ref. |fj| with the conclusion that 
it is consistent with all extant experimental and obser- 
vational constraints, provided > 0.24 eV. However, 
despite providing a natural explanation of the smallness 
of neutrino masses, this approach sheds no light on the 
observed pattern of neutrino mixing, nor on the hierarchy 
- if any - between the neutrino masses. These supposedly 
arise from the unknown structure of the Yukawa coupling 
matrix h in generation space. Any observed pattern of 
mixing and mass hierarchy among the neutrinos can be 
accommodated by this unknown structure, with the only 
constraint arising from the strict absence of neutrinoless 
double /3-decay. A simple generalization of this picture 
with heavy right chiral Majorana neutrinos in the bulk 
leads to light physical Majorana neutrinos via the seesaw 
mechanism. There, one also needs to invoke (unknown) 
flavor dependent brane-bulk couplings. 

There are higher dimensional scenarios 0, IE 0- E3 
with a somewhat opposite approach. These assume tiny 
flavor diagonal Majorana masses for left chiral neutrinos 
on the brane without explaining their smallness, but in- 
voke extra dimensions to generate their mixing pattern. 
Here the higher dimensional KK modes induce a seesaw 
leading to light physical Majorana neutrinos with flavor 
mixing. However, this mixing pattern is strongly con- 
strained. Suppose [HI that (1) no flavor dependence is 
assumed in the couplings between the brane and bulk 
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states and (2) flavor mixing is not allowed on the brane 
so that the brane Majorana neutrinos are kept flavor di- 
agonal 0, IE EH • Then, our observation is that the 
Dienes-Sarcevic model 7. 12], with three flavor diagonal 
neutrinos on the brane, is ruled out by the data. The 
allowed parameter space actually gets quite constrained 
@ even when an arbitrary flavor dependence is included 
in the brane-bulk couplings. We also find that a version 
of the model, given in Ref. [j| with four flavor diago- 
nal brane Majorana neutrinos and which has flavor blind 
brane-bulk couplings, is excluded. 

One attractive feature of the D-S model pfl is the abil- 



ity to provide a theory of neutrino mixing instead of 
just ascribing it to an unknown flavor dependence of the 
brane-bulk couplings. The intrinsic flavor diagonal Ma- 
jorana masses on the brane for the active neutrinos Uf (/ 
=1,2,3) are taken to be mj, while the brane-bulk link- 
age is characterized by a universal mass m appearing in 
all brane-bulk off-diagonal elements of the full neutrino 
Majorana mass matrix M.. Compactification to four di- 
mensions yields the zero mode as well as to a tower of 
sterile neutrino KK states, as mentioned above. After 
the imposition of appropriate orbifold conditions, the full 
neutrino Majorana mass matrix in the D-S model reads 
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In consequence of Q), the active neutrinos on the brane 
mix both among themselves and with the sterile KK 
states. Thus they undergo flavor oscillations on prop- 
agation. Let U be the unitary matrix which acts on 
mass eigenstate neutrino fields to yield flavor eigenstate 
neutrino fields. U T A4U is then diagonal. The matrix U 
therefore is given by the inverse of the matrix of eigenvec- 
tors of M and is computable as such. The time-averaged 
probability f/->/' that \vf) oscillates into \ vf), with both 
f and /' being flavor indices on the brane, is given by 

N i-1 

PjZf = S ft . MUfiU* ri U* fj U rj ], (2) 

where i and j are neutrino mass eigenstate indices. 

Let us now discuss the experimental results for the 
physical neutrino mass squared differences and the aver- 
age neutrino oscillation probabilities. The atmospheric 
neutrino data yield [ljj 13 

1.3 x 1CT 3 eV 2 < Am 2 tm < 3.2 x 1CT 3 eV 2 , (3) 
0.9 < sin 2 20 atm < 1. (4) 

The combined solar neutrino observations plus the Kam- 
LAND data imply [IE EE E3 

6.1 x 1CT 5 eV 2 < Am 2 Q i < 9.0 x 10" 5 eV 2 , (5) 
0.71 < sin 2 2(9 so i < 0.91. (6) 



The data [T£| from the CHOOZ reactor require that [3 

sin 2 26» 1 3 < 0.26. (7) 

In an effectively two flavor oscillation with a mixing angle 
<f>, the time-averaged oscillation probability is isin 2 20. 
We can therefore use the estimates 

0.45 < i^:<0.5, (8) 
0.355 < ?U < 0.455, (9) 
< 0.13 (10) 

in our analysis [2(i|. Finally, the recent WMAP result 
puts [2l| an upper bound on the sum of the physical 
masses of stable light neutrinos that mix among one an- 
other, namely 

^TO; < 0.71 eV. (11) 

i 

A constraint, therefore, on any model is that there must 
be at least three light neutrinos with the sum of their 
masses not exceeding 0.71 eV. 

It has been argued by Pierce and Murayama j2^| that 
the WMAP upper bound should include any light sterile 
neutrino that mixes with the active ones. This mixing is, 
of course, restricted to be rather small by the LEP con- 
straint on invisible Z decay. A wide window for the ex- 
istence of a sterile neutrino is nonetheless allowed by the 
WMAP results. Suppose the active fraction of the sterile 
neutrino is e so that e 2 ~ sin 2 (28 m i X ) , yielding e < 10 -1 . 
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Such a neutrino would have decoupled in the early uni- 
verse at a temperature T where G 2 F T 5 e 2 ~ T 2 Mp[. This 
temperature simply must be less than that of the QCD 
phase transition in order to avoid the dilution of those 
neutrinos by the entropy produced at this transition. 
Thus e must lie in the the range 1CP 6 < e < 10 _1 . This 
means that in most regions of the parameter space of the 
models examined here, the summation in includes 
the contribution of light sterile neutrinos. 

There are five parameters in the Dienes-Sarcevic (D- 
S) model : the input Majorana masses 7711,2,3 for the 
three neutrinos on the brane, the universal brane-bulk 
off-diagonal coupling m and the radius of the compact 
extra dimension R. For the present study we can set 
all mass scales in eV and R in eV" 1 . Let us also com- 
ment on the infinite number of sterile neutrinos in the 
model. There are, of course, certain limiting regions of 
the D-S parameter space when all the sterile neutrinos 
are sufficiently heavy to be effectively decoupled from 
the phenomena leading to © - . There are other re- 
gions where one or more sterile states can become part of 
the light neutrino spectrum. The existence of such light 
sterile neutrinos is already somewhat disfavored [23| by 
analyses of existing neutrino data. Though the situation 
is far from definite, results from the ongoing mini-Boone 
experiment are expected to provide a clearer conclusion. 
At the moment, several light sterile neutrinos can still be 
accommodated. 

We now come to a numerical investigation of the D-S 
model in the five-dimensional parameter space (mi, 777,2, 
777.3, m i R^ 1 )- An initial requirement is that there must 
be at least three light neutrinos with the sum of their 
masses obeying Thus those regions of the parame- 

ter space, which are not consistent with this constraint, 
are immediately excluded. We attribute the masses /zi.,2,3 
to the physical light neutrinos, always ordered [24[ such 
that 113 > /.12 > Ml- The conditions J2J and JSJ) then 
apply on the pertinent differences of squared mass eigen- 
values. After employing these conditions, we have made 
a complete scan of the parameter space. Our quest is 
to determine if there is any region in which all the other 
constraints, i.e. J7J) to iJTTjl. can be satisfied. Our answer 
is in the negative, i.e. there is no point in the parameter 
space where all the four conditions can hold simultane- 
ously. As will be shown below, there are parts of the 
parameter space where (JHJ), <|10|) and can be simulta- 
neously satisfied. Similarly, there are other regions where 
the conditions of ©, (fTUf> and (fTT|> can be met. But these 
two sets of points do not overlap. Our conclusion is that 
the D-S model is ruled out by the data. 

From practical considerations in making our scan, we 
have taken the number of contributing KK states N in @ 
to be finite though there actually are an infinite number 
of states in the KK tower. We choose N = 40 but have 
checked that varying this truncation affects our results 
negligibly. Furthermore, the choice of parameters in our 
numerical scanning was dictated by the eigenvalues of the 
mass-matrix (QJ, as described below. Consider, for the 



moment, the possibility of flavor dependent couplings for 
the brane neutrinos, leading to off-diagonal elements m\ 
of Ai. The eigenvalues of M. are then determined by the 
characteristic equation 

1 3 ™'2 E>2 

-tan ttA = J -. 12 

i—1 

Here, A = /iR and \x represents the eigenvalues. Let us 
also define d 2 = R 2 {J2i=i m 'i 2 ) so tnat m 'iR = dei, where 
J2 i=1 e| = 1. The two zeros of 

i— 1 

are given in the limit m\ — m! 2 — m' 3 — m (i.e., e\ = 

e2 = e 3 = e) by 

p = 2e 2 (mi + m 2 + m 3 ), 

1 /2 

C = [p — 4e (77727773 -I" 7773mi + 77117772)] . (14) 

It has been shown in Ref. that, as d increases, two 
of the eigenvalues move towards fib and fi c whereas the 
other eigenvalues move towards half-integral values times 
R . For tl > 1, the two eigenvalues sit at /if, and fi c 
(independent of R) , while the rest are fixed at those half- 
integer values times We have verified this property 
numerically. We have thus not included large values of 
to (mR = de) in our numerical scanning of the param- 
eter space, since that would not allow for two different 
mass squared differences as required to explain both the 
solar and the atmospheric neutrino deficits. The other 
significant constraint on the range of 7774,2,3 is tU - Let 
us now consider variations in the radius of compactifi- 
cation R. One should note first that, so long as the di- 
mensionless products miR, mi R, m^R and mR remain 
the same, one would get the same time-averaged oscil- 
lation probabilities. We have further checked that, for 
a fixed set of (7771,7772,7773 and m), our conclusions do 
not change with the variation in R permitted by the con- 
straints given in © and (|llfl . This way R can be 
seen to be constrained in the range aproximately given 
by 1(T 6 eV -1 < R < 10 eV -1 . 

In our convention [24[ , Am 2 ol = fi 2 , — fi 2 and A777^ tm = 
A*i — M2 ; f° r a normal ordering of neutrino masses, while 
^ m soi = A*3 ~ an d A TO atm = fj, 2 — fi\ in the case of 
an inverted ordering. Allowed regions, in the m\ — 7773 
plane for the choice of = 1 eV, are shown in Figs. 
(a) and (b). The incompatibility of the D-S model with 
the presently known facts from neutrino oscillation and 
cosmological data is made evident by these figures. In 
both figures mi, 777,2 and 7713 are taken to vary between 
0.001 eV and 0.31 eV, whereas m is varied between 0.001 
eV and 0.075 eV. Though, in the inverted few 
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FIG. 1: Allowed regions in the mi and m.3 plane for (a) the normal ordering and (b) the inverted ordering of neutrino masses. 
Choices of parameters are described in the text. 
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FIG. 2: Time-averaged oscillation probability (a) P a tm (in 
Blue/star) and (b) P so i (in orange/triangle) for the inverted 
ordering of neutrino masses as a function of 777,3. Only 
those points are considered here which satisfy for (a) the So- 
lar+CHOOZ and for (b) Atmospheric + CHOOZ constraints. 



points satisfying the conditions of @, i|l(J|) do lie close 
to the points satisfying (jHJ, l]lUp. in fact there is not a 
single point in the five-dimensional parameter space of 



(mi, m2, TO3, m, which is compatible with all the 

constraints. To further convince the reader, we display 
in Fig. [3 the time-averaged oscillation probability (a) 
Patm and (b) P so i as a function of 7713 in the inverted 
case for all points shown in Fig^ Clearly, there is no 
consistent region of parameter space. We have checked 
that this holds as R is varied in the range given above. 

Next, we consider the model of Lam 9] which requires 
a sterile neutrino (with a Majorana mass 7714) on the 
brane in addition to three active neutrinos. In the limit of 
large coupling between the brane neutrinos and the bulk 
(i.e., d ^> rriiR ,1) one can have three "isolated" eigen- 
values which can give rise to two different Am 2 , whereas 
the other eigenvalues sit at half-integer values. Following 
Lam, we also arrange the Majorana masses of the brane 
neutrinos in such a way that 

< mi < Hi < mi < (i2 < < < 7774, (15) 

where fii are the isolated eigenvalues. Furthermore, we 
assume that the rriiR do not have integral or half-integral 
values. We take a flavor-blind brane-bulk coupling m as 
in the earlier scenario. After a thourough scan in the 
parameter space with the above mentioned constraint in 
()15JI . we find that although in some part of the parameter 
space it is possible to satisfy the constraints in @ and 
(JSJ), there is no region where both the constraints in JSJ) 
and © can be satisfied either with a normal or with an 
inverted ordering of the neutrino masses. 

In summary, we are able to exclude the Dienes-Sarcevic 
model and that version of Lam's model which has 
flavor blind off-diagonal mass terms, by using all neutrino 
oscillation data and the WMAP constraint on the sum of 
neutrino masses. Extra dimensional models of neutrino 
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masses and mixing angles need flavor dependent brane- 
bulk couplings. 

P.R. and S.R. acknowledge the hospitality of the The- 
ory Group at SLAC where this work was initiated. We 
thank S. Choubey, S. Goswami, A.S. Joshipura and H- 



J. He for their helpful comments. The research of S.R. 
was supported by the Academy of Finland (Project num- 
ber 48787). S.R. would also like to thank the Lady Davis 
Fellowship Trust in Technion, Israel, for financial support 
during the initial stages of this work. 



[1] N. Arkani-Hamed, S. Dimopoulos and G.R. Dvali, Phys. 
Lett. B429, 263 (1998); I. Antoniadis, N. Arkani-Hamed, 
S. Dimopoulos and G.R. Dvali, ibid. B436, 257 (1998); 
L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 

(1999) ; ibid. 4690 (1999). 

[2] K.R. Dienes, E. Dudas and T. Gherghetta, Nucl. Phys. 
B557, 25 (1999); G.R. Dvali and A.Y. Smirnov, Nucl. 
Phys. B563, 63 (1999); A.E. Faraggi and M. Pospelov, 
Phys. Lett. B458, 237 (1999), J. Maalampi, V. Sipi- 
lainen and I. Vilja, Phys. Lett. B512, 91 (2001); Phys. 
Rev. D67, 113005 (2003); R.N. Mohapatra and A. Perez- 
Lorenzana, Nucl. Phys. B576, 466 (2000); ibid. B593, 
451 (2001); R. Barbieri, P. Creminelli and A. Strumia, 
Nucl. Phys. B585, 28 (2000); A. Lukas, P. Ramond, 
A. Romanino and G.G. Ross, Phys. Lett. B498, 230 
(2001), D.O. Caldwell, R.N. Mohapatra, and S.J. Yellin, 
Phys. Rev. D64, 073001 (2001); N. Cosme, J.-M. Frere, 
Y. Gouverneur, F.-S. Ling, D. Monderen and V. Van 
Elewyck, Phys. Rev. D63, 113008 (2001); K. Agashe and 
G.H. Wu, Phys. Lett. B498, 230 (2001); A. Ioannisian 
and A. Pilaftsis, Phys. Rev. D62, 066001 (2000); G. 
Bhattacharyya, H.V. Klapdor-Kleingrothaus, H. Pas and 

A. Pilaftsis, Phys. Rev. D67, 113001 (2003); A. Dighe 
and A.S. Joshipura, Phys. Rev. D64, 073012 (2001), N. 
Arkani-Hamed, S. Dimopoulos, G.R. Dvali and J. March- 
Russel, Phys. Rev. D65, 024032 (2002); T. Appelquist, 

B. Dobrescu, E. Ponton and H. Yee, Phys. Rev. D65, 
105019 (2002). 

[3] There are more complicated models consisting of right 
chiral neutrinos on the brane and / or left chiral ones in the 
bulk with warped geometric configurations 4|, but the 
use of flavor dependent brane-bulk couplings to explain 
the observed pattern of neutrino mixing is quite generic 
in these. 

[4] Y. Grossman and M. Neubert, Phys. Lett. B474, 361 

(2000) ; G. C. McLaughlin and J.N. Ng, Phys. Lett. 
B493, 88 (2000); S.J. Huber and Q. Shafi, J. Phys. G29, 
1807 (2003) and Phys. Lett. B583, 293 (2004). 

[5] S.M. Bilenky and S.T. Petkov, Rev. Mod. Phys. 59, 671 
(1989). 

[6] H. Davoudiasl, P. Langacker and M. Perelstein, Phys. 

Rev. D65, 105015 (2002). 
[7] K.R. Dienes and I. Sarcevic, Phys. Lett. B500, 133 

(2001) . 

[8] C.S. Lam and J.N. Ng, Phys. Rev. D64, 113006 (2001). 

[9] C.S. Lam, Phys. Rev. D65, 053009 (2002). 
[10] Q-H. Cao, S. Gopalakrishna and CP. Yuan, 

hep-ph/0312339 
[11] Thus the entire flavor mixing among the light physical 



neutrinos is induced by the KK modes. 

[12] An extension of the model of Ref. with three different 
off-diagonal mass elements, was investigated in the limit 
of their very large values by Lam and Ng and shown 
to be incompatible with the data. Our study, however, 
covers the entire parameter space of the Dienes-Sarcevic 
model with a uniform off-diagonal mass. 

[13] We use 90% C.L. limits throughout. 

[14] M. Shiozawa, super-Kamiokande coll. in XX Intl. Conf. 
Neutrino Physics and Astrophysics, Munich, May 2002 
(http://neutrino2002.ph.tum.del; S. Fukuda et al. super- 
Kamiokande coll. Phys. Rev. Lett. 81, 1562 (1998); ibid. 
85, 3999 (2000). 

[15] Q.R. Ahmad et al, SNO Coll., Phys. Rev. Lett. 87, 
071301 (2001); ibid. 89, 011301 (2002); S. Fukuda et al, 
super-Kamiokande coll. Phys. Lett. B539, 179 (2002); 
B.T. Cleveland et al., Astrophys J. 496, 505 (1998); W. 
Hampelet et al. Phys. Lett. B447, 127 (1999); J.N. Ab- 
duradutov et al. JETP 95, 181 (2002); Y. Fukuda et al, 
Phys. Rev. Lett. 77, 1683 (1996); M.H. Ahn et al, Phys. 
Rev. Lett. 90, 041801 (2003 ). 

[16] S.N. Ahmed et al, nucl-ex/0309004 

[17] M.C. Gonzalez-Garcia and C. Pena-Garay, Phys. Rev. 
D68, 093003 (2003). 

[18] M. Apollonio et al., Phys. Lett. B466, 415 (1999). 

[19] M. Maltoni, T. Schwetz, M.A. Tortola and J. W.F. Valle, 
Phys. Rev. D68, 113010 (2003); G.L. Fogli, E. Lisi, A. 
Marrone, D. Montanind, A. Palazzo and A.M. Rotunno, 
Phys. Rev. D69, 017301 (2004); A. Bandyopadhyay, S. 
Choubey, S. Goswami, S.T. Petkov and D.P. Roy, Phys. 
Lett. B583, 134 (2004). 

[20] Our conclusions follow just from these averaged proba- 
bilities without any need to use the detailed shape of the 
corresponding distributions. 

[21] D.N. Spergel et al, Astrophys. J. Suppl. 148, 175(2003). 

[22] A. Pierce and H. Murayama, Phys. Lett. B581, 218 
(2004); G. Bhattacharyya, H. Pas, L. Song, and T. 
Weiler, Phys. Lett. B564, 175 (200 3); P. Crotty, J. Les- 
gourgues and S. Pastor, hep-ph/0402049 

[23] M. Maltoni, T. Schwetz, M.A. Tortola and J.W.F. 
Valle, Nucl. Phys. Proc. Suppl. 114, 203 (2003); M.C. 
Gonzalez-Garcia, M. Maltoni and C. Pena-Garay, Phys. 
Rev. D64, 093001 (2001). 

[24] This convention of always keeping fj,3 > fj,2 > fii, which 
is different from the more common PMNS convention 
of attributing /Lti,2,3 to |^ e , M , r ) in the 'no-mixing' limit 
U — > /, has been followed for search convenience and is 
inessential to the derivation of the final conclusion. 



